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THE INVESTIGATION OF ADVANCED REMOTE SENSING TECHNIQUES FOR THE 
MEASUREMENT OF AEROSOL CHARACTERISTICS 

By 

Adarsh Deepak^ 


INTRODUCTION 

This is the final report on investigations performed under research 
grant NSG 1252. The overall aims of the research program proposed under 
this grant were to conduct investigations of advanced remote sensing tech- 
niques and inversion methods for the measurement of characteristics of 
aerosol and gaseous species in the atmosphere. Of particular interest 
were the physical and chemical properties of aerosols, such as their size 
distribution, number concentration, and complex refractive index, and the 
vertical distribution of these properties on a local as well as global scale. 
In this connection, the principal investigator. Dr. Adarsh Deepak, provided 
the necessary science support to the NASA/Langley Research Center (NASA/LaRC) 
technical monitor on numerous projects outlined below, which involved develop- 
ment -of remote sensing techniques for monitoring of tropospheric aerosols as 
well as satellite monitoring of upper tropospheric and stratospheric aerosols, 
development of computer programs for solving multiple scattering and radia- 
tive transfer problems, as well as inversion/retrieval problems. A necessary 
aspect of these efforts was to develop models of aerosol properties. 

The research investigations have been categorized under the following 
topics, which are. listed in the order in which they will be discussed in 
this report: 

I. The Solar Aureole Method for Determining Aerosol Characteristics 

II. Participation in University of Arizona- Aerosol and Radiation 
Experiment (UA-ARE) , May 1977 


^ Research Associate, Old Dominion University Research Foundation, P. 0, Box 
6369, Norfolk, Virginia 23508, and Research Associate Professor, Department 
of Physics and Geophysical Sciences, Old Dominion University, Norfolk, 
Virginia 23508. 




III. Multispectral Solar Radiometry for Determining Aerosol Size 
Distributions 

IV. Aerosol Physical Properties of the Stratosphere (APPS) Experiment 
(The Draper Laboratory AAFE Project) 

V. Stratospheric Aerosol Photographic Experiment (SAPE) 

VI. The UCLA Tropospheric Polarimeter and Multispectral Radiometer Project 
(The Polarimeter and the Multispectral Radiometer as Remote Probes 
of Aerosols) 

VII. International Interactive Workshop on Inversion Methods in Atmospheric 
Remote Sounding, Williamsburg, VA, December 15-17, 1976 

VIII. Modeling of Aerosol Properties 

This final report briefly describes some of the salient features of the 
.research investigations conducted under grant NSG 1252 for the period 
September 1975 to December 1977 and, in addition, identifies some important 
problems in which further research should be continued in the near future. 
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I. 


THE SOLAR AUREOLE METHOD FOR DETERMINING AEROSOL CHARACTERISTICS 


INTRODUCTION 

The solar aureole method has been successfully used in recent years to 
determine size distribution and concentration of atmospheric aerosols, particu- 
larly tropospheric aerosols. A long-range objective of the research performed 
in this method was eventually to develop satellite- or Shuttle-based remote 
sensing systems for monitoring atmospheric aerosols. 

The solar aureole is the area of enhanced brightness closely surrounding 
the Sun's disk (within about 20 degrees) caused mostly by aerosol scattering 
of sunlight. Since the aerosols scatter predominantly in the forward direction, 
the contribution of atmospheric haze to sky radiance for angles close to the 
Sun is roughly 10^ to 10^ times the contribution by molecular scattering; 
this is illustrated in figure 1. It is to take advantage of this large 
signal range that a simple, portable, photographic solar aureole measurement 
(PSAM) technique was developed at the University of Florida in 1970 (ref. 1) 
and has since been used to diagnose the aerosol size-altitude distributions 
(refs. 2 to 5) , by using the aureole radiance measurements along the almucantar. 
The almucantar is a conical scan of constant solar zenith angle with the local 
zenith as the axis. A solar aureole measurement program was subsequently initi- 
ated at NASA/Langley Research Center in 1974, and the photographic solar aureole 
isophote (PSAI) method (refs. 6 and 7) was developed. Isophotes are lines or 
curves of equal radiance. The PSAI method is an extension of the solar 
almucantar aureole radiance (SAAR) method. In the latter, the radiance 
measurements taken along the almucantar are used to infer the aerosol prop- 
erties. In the PSAI method, in addition to making the almucantar measurements, 
one takes advantage of the fact, which emerged from our computer studies, that 
the shape of the solar aureole isophotes is sensitively dependent on the 
characteristics of the aerosol size distribution. Suggestions for the use 
of solar aureole measurements to determine aerosol properties had also been 
made earlier by Deirmendjian (refs. 8 and 9) and other researchers (refs. 10 
to 14). Solar aureole measurements, taken with scanning photometers for 
determining aerosol properties, have also been made by Shaw (ref. 15) and 
Twitty et al. (ref. 16). 
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The following sections briefly describe the theory, photographic measure- 
ments, and results of the PSAI method. For details of some of the theoretical 
aspects discussed here, see references 1 and 17. 

SIMPLIFYING ASSUMPTIONS 

The forward- scattered radiance, being highly sensitive to n(r), is 
relatively insensitive to effects due to aerosol refractive index, polariza- 
tion, and multiple scattering (MS), a fact that helps in simplifying the 
inversion problem. In this paper, only the single scattering (SS) theory 
treatment is considered, which should help in understanding the difficulties 
involved in the inversion of aerosol scattering measurements. Therefore, 
one makes the following reasonable simplifying assumptions: 

1. Par.ticles are spherical; hence Mie theory results can be used 
in computations. 

2. ^ The atmosphere is horizontally homogeneous and vertically 

inhomogeneous . 

3. Absorption effects are ignored by selecting to work in spectral 
regions for which atmospheric absorption is nil. 

4. The polarization effects are small for forward -scattered light and 
.can be ignored. 

5. For relatively clear days (visibility greater than 15 km), the MS 
effects at the forward-scattering angle are small compared with SS 
effects (ref. 4) and can be ignored. 

6. An average value for the refractive index of all atmospheric 
aerosols is assumed for forward scattering. 

7. The atmosphere is treated as plane-parallel; the spherical Earth 
effects, which become significant for zenith angles (j> larger than 
75 degrees, are incorporated into the theory by the use of the 
generalized Chapman- type functions S((j)) (refs. 4 and 18) in place 
of the secant functions. 

The case of retrieval of size distribution from solar aureole data in 
the presence of multiple scattering has been discussed as part of another 
contract report (ref. 19) and, therefore, will not be discussed here. 
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SINGLE SCATTERING THEORY OF THE SOLAR AUREOLE 


Figure 2 illustrates the geometry of the calculation. Shown are an 
acceptance cone do originating at the detector and a solid angle dO' 
centered at the elemental scattering volume dV at altitude y(km); (|> 
and are the zenith angles of the Sun and the narrow view cone, and oi 
is the dihedral angle between the normals to the Sun zenith and view cone 
zenith planes intersecting at dV. The scattering angle \p is then given 
by the relation 

cos Tp = cos 4) cos <|) sin 4> sin 4> cos w (1) 

s s 

The element dV is given by 

dV = R^ dOS((j)) dy (2) 

where the generalized Chapman- type functions (refs. 4 and 17) 


S((j)) = sec <|> • (for ^ 75 degrees) (3) 

We make measurements in regions of the sky for which 4i < 75 degrees and 
_< 75 degrees; then, dV = R^dO sec 4> dy. 

The optical depth defined by 


Tj (X,y) 


=/' 


(x,y) 


(j = M,P) 


(4) 


for a ray traversing the distance from the Sun to the air mass element dV 
is given by 


Tj = S^Tj(X,y) sec (j = M,P) 


(5) 


and from the air mass to the detector by 
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■^2 = 2jjT^(X,0) - T^(X,y)j sec tj) 


( 6 ) 


where M denotes air molecules; P, the particulates (aerosols); X, the 

wavelength; and 3' the volume scattering coefficient (VSC) (km“^) at altitude 

y for the jth constituent. In this paper, all quantities represented by 

Tpj 0p> Sp, Fp, Fp, Pp, and Pp are functions of m (= m* - im") , 

the aerosol complex refractive index. The primes denote the y-dependence 

of the quantities. The terms 3* and 3,*, are defined as 

p M 

3p(X,y) =/ n(r,y)'irr2Q(x,m)dr _ ( 7 ) 

*'^2 


where Q(x,m) is the efficiency factor (ref. 20), x = 2iTr/X is the particle 
size parameter, r]_ and r 2 are the minimum and maximum ' values of r, 
and 

3j!iCX,y) = gj^(X)Pj^(y) (8a) 

where the VSC for molecules is 




8TT^(n^ - 1)^ 
NX'^ 


(4 + 3d) 
(4 - 3d) 


(8b) 


In equation (8b), N is the number of molecules (cm"^); n, the refractive 
index of the medium, d = 4A/ (1 - A) ; and A is the depolarization of scattered 
light at a scattering angle of 90 degrees for a linearly polarized incident 
radiation with its electric vector perpendicular to the scattering plane. 

For unpolarized incident light, A is replaced by ? = 2A/(1 + A). Then, the 
volume scattering function (VSF) for air molecules (j = M) is 

Fj^(i|j,X,y) = 3^j(X)Pj^(t)Pj^(y), km'^sr-l (9) 


where the molecular phase function is 

^ li? + cos2 ,|j), sr-1 


( 10 ) 


6 


where Pj,|Cy) is dimensionless function representing the altitude dis- 
tribution of molecular density. The VSF for aerosols (j = P) is 


Fp(ijj,X,y) = 3pCX,y)Pp(iJ<,X,y) 


( 11 ) 


where the aerosol volume phase function is 

= -^ — n(r,y) jii(i(j,m,x) + i 2 (i|',m,x) j dr (12) 

P 2k2B;(X,y3 4i ' ^ 

I 

and ii and i 2 are the Mie intensity functions and k = 2 tt/X. 

The sky radiance due to the molecules and aerosols in the volume element 
dV is then given by 

dB((j),(j)^,co,X) = H^(X) sec ^ (FjJj(ij),X,y) 

, . -(T1+T2) 

+ Fp(ii),X,y)j e dy (13) 

Integrating over all such elemental volumes, the total single scattered sky 
radiance is 

I " 

B(4>,<i>g,to.X) = gJ^ {Fj^(i|;,X,y) + F’(i!j,X,y)j dy (14) 

where 


G = H^(^) sec ^ e 


-x(X,0)sec (j) 


and 


D = sec - sec (}> 


(15) 


(16) 


Before discussing the inversion problem, the photographic solar aureole 
measurement (PSAM) technique will be described. 
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V 


THE PHOTOGRAPHIC SOLAR AUREOLE MEASUREMENT TECHNIQUE 

A photographic technique of making measurements of the solar aurcoU 
radiance is briefly described here. Figure 3 schematically illustrates the 
equipment. Photographs of the Sun’s aureole are taken with a small format 
camera (35 mm or 70 mm) through a wavelength filter with the Sun occulted 
by a neutral density (ND) disc (ND = 4) held coaxially on a stem about 0.6 m 
to 1.3 m in front of the lens. The ND filter attenuates the radiance of 
the direct sunlight. .by a factor of lO'^, and the optical densities of the 
surrounding aureole are shown in a typical photograph in figure 4a. The Sun's 
image not only enables one to calibrate the entire photograph relative to its 
radiance, but also enables accurate measurements to be made of the angular 
distances from the Sun. Photographs are taken through different wavelength 
filters. In addition, the direct solar irradiance measurements through the 
same filters are made with a photometer to determine t(X,0).. 

A. Solar Aureole Isophotes 

I 

The optical density of the whole photograph is read with the help of a 
Joyce-Loebl Isodensitracer, which gives digital data output on.a magnetic tape 
and, at the same time, provides an isodensity tracing, such as shown in 
figure 4b. Isodensities are lines or curves of equal optical density in a 
photograph. Isophotes are then generated from the taped data output, as 
shown in figure 4c, where an economy-wise reduction has been made in the 

number of data points . 


B. Almucantar Radiance 

The photogrmmetry of the solar aureole is disoussed in detail in refer- 
ence 21. It is shown that the almucantar projects on the film as a conic 
Cfig. 5). The shapes of the conics for different values of the solar zenith 
angle 4> are illustrated in figure 6. Accordingly, the measured almucantar 
radiance^as a function of the scattering angle is shown in figure 7. The 
peak at zero degrees corresponds to the Sun's direct light reduced in inten- 
sity by the ND disc. With equation (14j , it is possible to obtain the infor- 
mation about the size-altitude distribution of aerosols from the multispectral 

measurements of the sky radiance, B(X), as a function of ♦, and m, 

and of total optical depth t(X,0). 
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Many techniques for obtaining the aerosol characteristics suggest them- 
selves in the light of this analysis. One of the simplest and most elegant 
is the method based on the almucantar radiance measurements, as explained in 
the following section. 

ALMUCANTAR RADIANCE AS AN INDICATOR OF SIZE DISTRIBUTION 

In the almucantar, which is a conical scan of constant solar zenith 
angle (j)^ with the local zenith as axis, the total radiance B is equal 
to the sum of radiance due to air molecules (Bj^) and particulates (Bp) , so 
that 

BpCi|>,X) = B - Bj^ (17) 

where ^ is the scattering angle, and X is the wavelength of solar radia- 
tion. From equation (17), one can then experimentally obtain the volume 
scattering function Fp(\|;,X) as a function of ip, the quantity Fp being 
a functional of n(r). It is the experimental Fp(4i,X) vs. ip data which 
is inverted to obtain n(r). The experimental data for the altitude-integrated 
VSF is then given by 

B(<j) = <|) ,0J,X) 

^ Pj^(4<,X) j^(X,0) (18) 

RETRIEVAL OF AEROSOL SIZE DISTRIBUTION FROM ALMUCANTAR MEASUREMENTS 

Retrieval aerosol size distribution from equation (17) can be accomplished 
by either the model-fitting approach or the numerical inversion scheme. 

A. Model Catalog Search Method 

In this approach, an analytic model for n(r) (ref. 22) is chosen from 
a parametrized catalog of phase functions which is generated for different 
values of m (ref. 23). As an illustration, the average value of m is 
assumed to be 1.55 - i(O.O), anci a typical model for n(r) is chosen as 
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(Pl/r)'’"* 

- 1 

where pi, P2, P3, and v are adjustable constants. Typical size distri- 

bution curves for different values of the parameter v are shown in figure 8. 
Any other analytic model of n(r) could do just as well. 

A set of phase function curves corresponding to values of v that range 
from 4.0 to 5.2 are shown in figure 9. By comparing the experimental phase 
function curve in figure 8 with one of the curves in the catalog, one obtains 
a reasonably good estimate of the n(r). If, however, one wants to go a step 
further and obtain the n(r) that gives the best fit to the experimental phase 
function curve, a nonlinear least squares computer code is used. This is 
described next. 



P3 

n(r) = -5- 


P2 


(. 


B. Numerical Inversion Scheme 

The experimental Fp(i|^,X) vs. tj; data was inverted by the use of non- 
linear least squares (NLLS) method. On the basis of the exact Mie theory, 
Fp(\|;,X) is defined by 


FpC<]^,Xj 




(il + i2) nCrjdr 


( 20 ) 


where ii and i2 are Mie intensity coefficients (ref. 18) r^ = 0.04 ym and 
r2 = X3 pm. The values i^ and i2 are computed using J.V. Dave's computer 
code. A two-term model of n(r) was selected with adjustable parameters p^, 
i = 1, 2, 3, and 4. Each term is a Haze M model (ref. 23), so that 


n(r) = pir {exp(-p2 v^) + P3 exp(-p4 v^)} 


( 21 ) 


Initial estimates of the parameters are obtained with the help of a parametrized 
graphical catalog of phase function plots (ref. 24) corresponding to different 
size distribution models (ref. 22). The best fit values of p^^^ are obtained 
by the use of the NLLS computer program. The method used in this program is 
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based on Taylor's series expansion o^ the function in terms of its param- 
eters about the current estimates of the parameters and iteration of the 
resulting linear approximation until the estimates converge in the least 
squares sense. Results of typical inversion of experimental data for X 
of 600 nm, 500 nm, and 400 nm, based on SS, are shown in figure 10. The 
experimental measurements were obtained during a field experiment described 
in section II. 

RETRIEVAL FROM SOLAR AUREOLE ISOPHOTES 

An extensive parametric computation of the solar aureole isophotes as 
functions of n(r) and Ppf.y)> 3-nd to some extent of m, has been carried 
out by using equations (14) and (19). Everything else being the same, circum- 
solar isophotes corresponding to three values of the parameter v of the size 
distribution n(r) (fig. 7) are shown in figures 11a, 11b, and lie. It is 
easy to see from figures 7 and 11 that, by slightly increasing the number of 
smaller particles and decreasing the number of larger particles as v becomes 
larger, the shape of the isophotes undergoes a dramatic change. As the value 
of V increases upwards from a value of 4.0, the isophote pattern attains a 
shape increasingly similar to that of the experimentally obtained isophotes 
(fig. 4c), until for a value of v = 5.0 the computer-generated pattern (fig. 
11c) best resembles the latter. Thus, figure 11a illustrates the fine sensi- 
tivity of the patterns of the circumsolar isophotes to the size distribution 

n(r) . 

■ Attempts are underway to optimize computer programs that will obtain the 
best fit to the isophote pattern or numerically invert the isophote data by 
utilizing, to their advantage, the fact that radiance along each isophote 
remains constant. In this regard, it is important to keep in mind that, 
even though isophotes depend on both n(r) and pp(y), within the aureole 
region the isophote shape is more sensitively dependent on n(r) than on 
p (y) . The sensitivity of the isophotes to Pp(y) increases for larger 
angular distances from the Sun. 

CONCLUDING REMARKS 

In contrast to the modeling approach, in which one starts with a model 
judiciously hypothesized from experience, one begins the numerical inversion 


approach with an initial estimate of parameters. A typical iterative scheme 
is described by Malchow and \Vhitney in reference 25. Simulations of sky 
radiance isophote shape or the phase function are performed by using equation 
(14) or equation (20), respectively, along with the initial set of parameters. 
Comparisons are made with the appropriate measurements. If the discrepancy 
is greater than a minimum prescribed value, then a parameter updating algorithm 
is used to obtain a new estimate of parameters. However, when the convergence 
criterion is satisfied, the final estimates of parameter are assumed to be 
accurate. This is very similar to the NILS scheme. 

It should be clearly pointed out that, since the radiative transfer 
program is called several times during each iteration in such an interative 
inversion scheme, use of the Mie theory computations of scattered radiance 
in an inversion scheme can often become prohibitively expensive. Thus, in 
developing a radiative transfer code for use as an inversion scheme, it is 
imperative that it be computationally as fast as possible in order to be of 
practical use in numerical retrievals of aerosol characteristics. 
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II. PARTICIPATION IN UNIVERSITY OF ARIZONA-AEROSOL AND RADIATION 

EXPERIMENT (UA-ARE) , MAY 1978 

The principal investigator, in concert with R.R. Adams of NASA/Langley 
Research Center, participated in a multi-university, multi -instrument coor- 
dinated aerosol and radiation experiment (ref. 26), conducted by Dr. J. Reagan 
at the University of Arizona in May 1977 . The solar aureole photographic 
equipment was taken to Tucson, Arizona, and over 500 photographs of the Sun's 
aureole were taken from the rooftop of the Engineering Building. Aureole 
photographs were also taken at various altitudes (2000 to 9000 ft; 609.6 to 
2,753.2 m) as we drove up to Mt. Lemon on the Catalina Mountains west of Tucson. 
In addition, radiometer readings of the direct solar radiation were taken through 
the same spectral filters -that were used for photography. 

These experiments were performed at the same time the other radiation 
measurements, such as Sun radiometers, diffuse radiation measurement instru- 
ments, bistatic lidar, were being deployed, thereby providing not only con- 
sistency checks but also intercalibration of the data. Aerosol and radiation 
measurements were also taken from aboard the Pennsylvania State University's 
instrumented plane. 

The nonlinear least-squares scheme described earlier was used to invert 
the solar aureole data for X = 400, 500, and 600 nm. The retrieved results 
are summarized as follows (ref. 19): 

The size distribution model was the sum of two Haze M models [eq. (21)]. 
Computations were performed in order to understand the effect of the following 
parameters, viz, the refractive index, integration limit r^ and, ignoring 
multiple scattering, the retrieved size distribution. The results are 
"described as follows: 

Three refractive indices were used, namely, 1.45 - i(O.O), 1.50 - i(O.O), 
and 1.50 - i(O.Ol). Essentially no difference was found. in the retrieved 
results obtained for the different refractive indices, although the param- 
eters were slightly different. In addition, the decrease of ri from 0.04 
ym to 0.01 ym made no difference in the final n(r) results. For scatter- 
ing angles greater than 12 degrees, it was found that multiple scattering 
should be taken into account to improve the accuracy of results. 
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Assuming single scattering, size distribution results were obtained 
by inverting the three sets of experimental data (fig. 10), first indi- 
vidually (ref. 26), and then simultaneously. Comparison among the retrieved 
size distributions obtained individually shows that the SD’s for X = 0.5 ym 
and 0.6 yra, are closer to each other than to SD for X'= 0.4. The main dif- 
ferences occur in the region around 1,0 ym, with SD for X = 0.4 ym having 
fewer particles (fig. 10). The correct method of obtaining the SD is by 
inverting the three sets of data simultaneously. The results thus obtained 
agree well with those obtained individually and also with the ground truth 
data obtained by airborne Whitby ((8» and Royco ( + ) counters, which have 
been plotted in figure 10. Even though the ground truth data is sparse 
and is for May 7, 1977, the day following May 6, 1977, the day on which 
solar aureole measurements were made, they were reasonably close to the 
retrieved results and gave us confidence in the latter. 


III. MULTISPECTRAL SOLAR RADIOMETRY FOR DETERMINING 
AEROSOL SIZE DISTRIBUTIONS 

Measurements of the total optical depth t(X) of the atmosphere were 
made by measuring the direct solar radiation at various wavelengths (X) by 
means of a solar radiometer. These measurements can be inverted — independently 
or in conjunction with solar aureole data — to retrieve aerosol size distribu- 
tion. We used a Tektronix radiometer (half-cone angle ~ 4 degrees) in our 
experiments. Whenever such a radiometer is used to measure the direct beam 
from a light source (the Sun in our case) , some diffuse sky radiation invari- 
ably enters the detector's field of view along with the direct beam. If 
is the unattenuated intensity and I(X) is the intensity after passing through 
a mediiam, then the optical depth x(X) is given by Bouguer's (or Beer-Lambert 's) 
transmission law 

I(X) = I^(X) 

which is valid only for the direct radiation. 

The experimentally measured intensity I (X) , however, contains not only 
the direct beam, but also some forward- scattered radiation. Therefore, the 
optical depth obtained by the application of Bouguer's law is strictly speaking 
not the true optical depth, but an apparent optical depth. Thus, a forward- 
scattering correction factor R is often used to account for this scattered 
radiation entering the detector's field of view. In order to understand the 
behavior of such a correction factor as a function of the particle size distri- - 
but ion and the complex refractive index, we carried out a theoretical parametric 
study first for the monodisperse aerosols (ref. 27) and then for polydispere 
aerosols (ref. 28). These two references are included as Appendixes 1 and 2 
of this report. The aerosol medium was assumed to be vertically and horizon- 
tally homogeneous. In this study, only single scattering was considered; multiple 

scattering was treated as being negligibly small (It should be pointed 

out here that the general case in which multiple scattering and vertically 
inhomogeneous atmosphere are present will be treated in reference 29 as 
part of another NASA contract.) 

Some of the conclusions of this investigation are summarized here as 
follows. It is shown that the forward- scattering correction factor, R, 
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decreases as the peak o£ the size distribution decreases. In addition, as 
the polydispersity o£ the size distribution becomes greater, the relative 
proportion o£ larger particles increases, resulting in an increase in the 
£orward scattering o£ radiation within narrower angles. Other conclusions 
are that the lens-pinhole type detector optics, in most practical cases, 
yields more accurate transmission/ extinction results than a detector without 
a lens-pinhole system. The conditions £or optimum designs o£ radiometer/ 
transmissometer detector systems are discussed in Appendixes 1 and 2 £or 
the cases o£ monodispere and polydisperse aersol media, respectively. As a 
typical example, i£ we de£ine error £actor E = (1 - R), then the percentage 
error contours representing trade-o££s between particle size parameter 
x(= 2irr/X) and the detector viewcone hal£-angle 6 degrees are shown in 
£igure 12 £or m = 1.50 - i(O.O). 
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IV. AEROSOL PHYSICAL PROPERTIES OF THE STRATOSPHERE (APPS) EXPERIMENT 
(THE DRAPER LABORATORY AAFE PROJECT) 

The principal investigator, Adarsh Deepak, assisted the technical monitor 
in Phases I and II of the APPS experiment developed under the AAFE project by the 
C.S. Draper Laboratory. Phase I involved the development of the theoretical 
algorithms and computer codes for inverting the spectral measurements, and 
Phase II involved the design of the instrument. 

The development of the APPS experiment was aimed at the determination of 
the vertical profiles of atmospheric constituents and aerosol physical prop- 
erties about 10 km altitude in the stratosphere and lower mesosphere. The 
experimental plan was to make simultaneous multispectral measurements of sun- 
light scattered by the Earth's atmospheric limb by a photometer at various 
narrow spectral bands from aboard the Space Shuttle platform, and to develop 
computer codes to invert these measurements for retrieving vertical profiles 
of the concentrations of Rayleigh scatterers, ozone, and aerosol extinction, 
aerosol size distribution, and real and imaginary parts of the complex refrac- 
tive index. (These results are described in reference 25 by Malchow and 
Whitney.) 

The theory for solving the radiative transfer equation for a spherical 
aerosol atmosphere and for the inversion problem is very complicated. However, 
theoretical studies performed by the author helped in redirecting some of the 
theoretical efforts being performed at the Draper Laboratory to yield satis- 
factory results within a brief period. The computer code, namely. Sunlit 
Limb- Inversion Code (SLIC), developed at Draper Lab, was. adapted for -use 
on the 'LaRG computers. The lengthy computer programs were first checked 
out and later modified to include plotting capabilities. Some of the bugs 
in the program were removed and the programs were made more sophisticated. 

A typical plot showing the inversion of the simulated data for 13 km to 
24 km for an eight-channel sunlit limb radiance profile to obtain aerosol 
extinction profile is shown in figure 13. In the figure, circles represent 
the true state; squares, the state based on initial guess values of the 
parameters; diamonds, the state after the first iteration; and triangles, 
the final state after 2 to 10. iterations at different altitudes. The 
convergence to the true state for a maximum of 10 iterations seems rather 
good. It should improve for a larger number of iterations. 
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V. STRATOSPHERIC AEROSOL PHOTOGRAPHIC EXPERIMENT (SAPE) 

The theory and inversion methods involved in SAPE are essentially the 
same as those in the APPS experiment, except that the multispectral measure- 
ments of the sunlit limb radiance profiles are obtained with a photographic 
camera rather than with a photoelectric radiometer. Whereas the A 
experiment Is based on the measurements through eight spectral channels 
located in the near UV, visible, and near IR. SAPE involves measurements 
of limb radiance and its polarisation through two wavelength filters. 

The overall objective of the SAPE is to demonstrate the capability of 
the method to determine the stratospheric aerosol characteristics Cphysica 
and spatial) from the measurement of the solar radiation scattered ^ 
Earth’s atmospheric limb when the Sun is in any known position using sophis - 
cated inversion techniques. If the Sun is within the limb, then t e ire. 
attenuated solar radiation is measured in addition to the forward-scatt 
radiation, so that inversion of the two sets of data yields results about 
aerosol characteristics that complement each other. It is the a sence o 
restrictions upon the Sun's location with respect to the observed portion 
of the sunlit limb that gives this technique a temporal and global coverage 
that is greater than that achieved by any other existing technique. 

The experimental objectives of SAPE are to measure, photographically, 
the profiles of (a) limb-scattered radiance and polarization, and (b) limb- 
attenuated direct solar irradiance and Sun- shape through narrow spec ra 
regions at 880 nm and 650 nm, when the Sun is at various spacecra 
zenith angles. The exposure times needed to take analyzable photographs 
of the sunlit limb through the two filters have been evaluated. The 
photogrammetry. sensitometry. densitometry, and the photometry needed to 
obtain accurate radiance data have been extensively investigated. 

computer programs for inversion schemes required for retrieving aerosol 
and gaseous characteristics have been optimized. The anticipated resu s 
SAPE are the profiles of stratospheric aerosol extinction, number density, 
ize distribution, co^lex refractive index, and atmospheric refractive index. 
In addition, it is anticipated that detection of high-altitude particula e 
Tayers at si and R5 km will be made possible. The one thing that we need 
now is to obtain one or more photographs of sunlit limb from aboard the 

Space Shuttle. 


18 


\ 


VI. THE UCLA TROPOSPHERIC POLARIMETER AND MULTISPECTRAL RADIOMETER PROJECT 
(THE POLARIMETER AND THE MULTISPECTRAL RADIOMETER AS 
REMOTE PROBES OF AEROSOLS) 

This project involved the determination of the tropospheric aerosol 
characteristics from the measurements of the polarization of scattered 
solar radiation in the Earth's turbid atmosphere. The measurements of 
scattered intensity and degree of polarization were performed with a 
NASA- funded polarimeter at UCLA by Dr. J. Kuriyan (ref. 30). It has been 
shown by several researchers that polarization of scattered radiation is 
sensitive to the aerosol refractive index. The inversion of the polarization 
data to obtain information about the aerosol properties is quite a formidable 
problem. So far, only the catalog search technique has been used (by Dr. J. 
Kuriyan) to retrieve aerosol SD and refractive index information. Numerical 
inversion schemes for polarization data have been investigated by the author, 
but a computer program has not yet been completed. The author plans 
to continue to develop such an inversion code. In addition, work on the 
sensitivity analysis of a Space Shuttle-based photopolarimeter technique 
was initiated at UCLA. The study was aimed at finding the optimum way in 
which the polarization measurements of upwelling radiation should be per- 
formed from aboard the Shuttle for global monitoring of the atmospheric 
aerosols. The algorithm for making such measurements and inverting the 
polarization data v^as based on the catalog search method. However, there 
is a strong need to develop a numerical inversion scheme for retrieving 
aerosol characteristics from polarization measurements of both upwelling 
and downwelling radiation. 
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VII. INTERNATIONAL INTERACTIVE WORKSHOP ON INVERSION METHODS IN 
ATMOSPHERIC REMOTE SOUNDING, WILLIAMSBURG, VIRGINIA, 

DECEMBER 15-17, 1976 

This workshop was conducted under contract/task authorization NASl-14193-24 
as a part of grant NSG 1252. It was organized and chaired by the principal 
investigator, Adarsh Deepak, in cooperation with LaRC Technical Monitor, 

Dr. M.P, McCormick, who was the Associate Chairperson. 

The ODU-Langley Research Center Workshop on Inversion Methods in 
Atmospheric Remote Sounding was a resounding success. Professor J, Lenoble, 
University de Lille, France, reported on the Workshop in Applied Optics 
(ref. 51). The workshop attendance, which was by invitation only, consisted 
of 73 participants from seven countries, representing universities, research 
laboratories, and U.S. Government agencies. The presentation of invited 
papers by some of the world's leading experts and the interactive discussions 
made the workshop a very stimulating and valuable experience for everyone. 

The workshop provided an interdisciplinary forum to review and assess 
the state-of-the-art in inversion methods available for retrieving information 
from remotely sensed data. The emphasis of the invited papers and the followup 
discussions covered the assumptions, methodology, resolution, stability, 
accuracy, and the future efforts needed in the various inversion methods. 

In this regard, the direct radiative transfer methods and results relevant 
to the inversion problem were also briefly reviewed in a special session 
on radiative transfer, held jointly with the OSA Topical Meeting on 
Atmospheric Aerosols, which preceded the workshop. This allowed researchers 
in different areas of remote sounding to compare and optimize the utilization 
of these inversion procedures in their respective remote sounding techniques. 
The workshop provided an avenue for free exchange of ideas and information 
among researchers concerned with the various inversion methods. 

Twenty-two invited speakers presented papers covering the mathematical 
theory of inversion methods and the application of these methods to the remote 
sounding of atmospheric temperature, relative humidity, and gaeseous and 
aerosol constituents. 

The invited papers and the recorded discussions from the radiative session 
and the inversion method sessions of the International Workshop were edited by 
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the principal investigator and published as workshop proceedings, entitled 
"Inversion Methods in Atmospheric Remote Sensing," by NASA and the Academic 
Press (ref. 32) . 
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VIII. MODELING OF AEROSOL PROPERTIES 

A parametric study to develop graphical catalogs of (a) various analytic 
methods for aerosol size distribution and cumulative size distribution and 
(b) aerosol optical properties was initiated under grant NSG 1252. Work on 
the first catalog has been completed (ref. 22) and will shortly appear as a 
NASA publication. It is intended to be of great help to aerosol researchers 
in representing their experimental data by commonly used analytic functions. 

The study further shows the equivalence of analytic functions in best fitting 
the experimental size distribution (SD) data. The work on the second catalog 
(ref. 23), viz, for the , optical properties associated with each of the size 
distributions is not yet complete. It will be completed under another contract. 
The computation of the optical properties is based on Mie theory computer codes. 
The optical properties of interest include the phase function and the extinction 
coefficient. 

Eight analytic models for n(r) were considered. Their properties are 
summarized in table 1. A typical parameterized graph of n(r) vs. r for 
model 2, viz. 


n(r) = ^ 


P3-1 

Pl\ C^/P2) 

P2/ 7 \P4 




(r/pz) 


P3, 


is shown in figure 14. 1"^ table 1, n(r) represents the size distribution 

(cm“^, um“^); N(r), the cumulative size distribution, cm"^; r^, the mode 

radius; M, , the kth moment; p. , i = l,2,...n the model parameters; and 
K 

the equation numbers mentioned in the table refer to the equations in 
reference 22 (to be published in the near future as a contractor report) . 


22 


CONCLUDING REMARKS 


Following are some of the important conclusions that emerged from the 
investigations described earlier. 

The solar aureole photographic method for determining atmospheric aerosol 
characteristics is a simple and practical method, whose validity has been 
well demonstrated. As a result of these investigations, a niunber of publica- 
tions have been issued by the principal investigator in co-authorshig with 
coworkers at the Institute of Atmospheric Optics and Remote Sensing and NASA/ 
Langley Research Center. In addition, several other groups around the world 
have started using this technique for the measurement of aerosol characteristics. 

Our results on the forward-scattering corrections to optimum extinction 
measurements have been successfully checked out by Professor Allan Carswell, 

York University. Toronto, Canada, and his graduate student (as mentioned to 
the author in a private communication) . These results were also helpful in 
interpreting NASA/LaRC's Tektronix radiometer data for solar radiation, 
obtained in conjunction with the solar aureole experiment. 

The SLIC code developed for inverting horizon radiance profiles to 
retrieve stratospheric aerosol and gaseous characteristics is an extremely 
fast code, which takes account of multiple scattering in a spherical atmos- 
phere. It needs further work done on it in order to make it more versatile 
so that extinction data as well as polarization data can also be inverted. 

Expertise has been gained in the various aspects of solar aureole 
photography, 'as well as in the* theory and inversion schemes investigated 
in connection with SAPE, to the point where now some photographs of the 
sunlit limb are needed from aboard the Space Shuttle. 

Polarimetry of the upwelling and downwelling radiation provides a 
powerful technique for determining the aerosol characteristics. In order 
to make it a practical remote sensing method, a fast computer code for 
numerically inverting polarization data needs to be developed. The principal 
investigator intends to continue to develop such a code, which is not too 
far from completion. 

There are three main components to remote sensing, namely, (a) the 
mathematics; (b) the physical concepts; and (c) the measurement techniques. 
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Much of the mathematics of the inversion problem were worked out in the 1960 's 
In order to make further improvements in accuracy of remote soundings, one 
needs to look more deeply into the physics and measurements aspects, 'ihe 
1./76 Williamsburg Workshop on Inversion Methods summarized many of the inver- 
sion techniques. It is suggested that a symposium or workshop dealing with 

the physics and measurement aspects of remote sensing be held in the near 
future. 

The parametrized catalog of aerosol size distributions has proved quite 
valuable in representing the experimental size distribution data by analytic 
models. The completion of the corresponding parameterized catalog of the 
optical properties will be very useful in various aspects of aerosol research. 
The principal investigator intends to continue to work toward its completion. 
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APPENDIX 1 


FORWARDSCATTERING CORRECTIONS FOR OPTICAL EXTINCTION 
MEASUREMENTS IN AEROSOL MEDIA. 1: MONODISPERSIONS 


Reprinted from APPLIED OPTICS, Vol. 17, page 2900, September 15, 1978 
Copyright 1978 by the Optical Society of America and reprint)^ by permission of the copyright owner. 


Forwardscattering corrections for optical extinction 
measurements in aerosol media. 1 : Monodispersions 


Adarsh Oeepak and Michael A. Box 


This paper, Part 1 of two papers, praaanu a parametric study of the forwardscattering corrections for experi- 
mentally measured optical extinction coefficients in homogeneous aerosol media, since some forwardacat- 
ured light invariably enurs. along with the direct beam, into the finite apetmre of the detector. Part 1 
treats the case of monodispersions: Part 2, that of poiydispersions. Forwardscattering is considered a sin- 
gio-scattering phenomenon, and the corrections are computed by two methods: one, using the exact Mis 
theory, and the other, the approxinute Rayleigh diffraction formula. A parametric study of the dependence 
of the corrections on the particle size parameter, real and imaginary parts of the complex refractive index, 
and the half-angle of the detector’s view cons has been carried out. The parameter ranges in which the re- 
sults obtained by the approximate formulation agree well with those obtained by the Mie theory are also in- 
vestigated. The agreement is especially good for small view cone angles and large particles and improves 
even more for slightly absorbing aerosol particles. Also discussed is the dependence of these corrections 
upon the experimental design of the transmission measurement systems. 


I. Introducilon 

The forwardscattering correction refers to the cor- 
rection that must be made to the optical extinction 
measurements carried out in a scattering medium, due 
to ths fact that along with the direct radiation lome 
forwardscattered radiation invariably enters^ihto the 
finite aperture of the detector. In fact, no known light 
measuring device exists that can measure the intensity 
of the direct beam of radiation, with the complete ex- 
clusion of the light scattered in a forward direction.^*^ 
In making extinction measurements, one may optimally 
design the experimental apparatus so that the effects 
due to the forwardscattered radiation (and stray light) 
are minimized,^ but still in order to get the true ex- 
tinction measurement, one must calculate and, if sig- 
nificant, correct for the forwardscattering contribu- 
tion. 

The term forwardscattering is considered a single, not 
multiple, scattering phenomenon in which scattered 
radiation reaches the detector after being scattered only 
once by scatterers situated within the path of the direct 
radiation. The case of multiple scattering, which causes 


additional corrections to the optical extinction mea- 
surements, will be discussed in a subsequent publica- 
tion. 

Two methods— one exact and one approximate — for 
deriving the forwardscattering corrections to the 
transmission law for both absorbing and nonabsorbing 
monodisperse aerosols are discussed in this paper. 
These methods have been discussed by Gumprecht and 
Sliepcevich^’^ for nonabsorbing monodisperse aerosols 
in connection with the problem of arriving at the opti- 
mum experimental design characteristics for an ex- 
tinction experiment. Excellent discussions of the for- 
wardscattering corrections are also given in Refs. 4-6. 
Here, however, the exact Mie thsory is given in a closed 
form solution, which increases not only th^ accuracy of 
the results but also the efficiency of the computations. 
The approximate method is based on the Rayleigh 
diffraction theory for large particles. Approximate 
expressions for the correction factor at the small Ray- 
leigh particle limit are also given. The discussions in 
this paper are confined to homogeneous spherical 
aerosols, having uniform size distribution, concentra- 
tion, and composition, such as may occur along hori- 
zontal paths in field or laboratory situations. 


The authors are with Institute for Atmospheric Optics & Remote 
Sensing, P.O. 3ox P, Hampton, Virginia 23666. 
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II. Transmission Law. (Bouguer's Law) 

An electromagnetic plane wave of wavelength X auid 
intensity 7o(X) (W/cm--Mm), after traversing a distance 
L through a homogeneous aerosol medium (Fig. 1), has 
intensity /(X), given by Bouguer’s Law 
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/(X) - /o(X) #xp(-r(X)l, (1) 

where r(X), the optical thickness, is defined as 

WX)-^(X)L, (2) 

and the total volume extinction coefficient d.,t(X) is 


Qtcat ia based on the total amount of light lost from the 
^am by scattering by a particle in all directions, and 
V»c«tf on the amoimt of light lost by scattering in all di- 
rections except within a cone of half-angle 6 in the for- 
ward direction, so that their difference is given by 


^•it(X) • fittuiW + [cm“‘], (3) 

the subscripts ext, scat, and abs denoting extinction, 
scattering, and absorption. For spherical monodisper- 
«ons of number density Np{r) (cm-3) and complex re- 
fractive mdex m » m' - im', 

« Tr*Qj{xjn)Np{r), j - ext, icat, and aba, (4) 

where x » 2rr/\ is the size parameter, and Q., the effi- 
cient factors. If the molecular contribution to d,., is 
negligible, from Eqs. (3) and (4) one obtains 

Qm • QtcM + Qite. (5) 

where 


■A f'«. 

X* Jo 


it) sinfidfi. 


( 6 ) 


and ii and ij are functions offf.x, and m, called the Mie 
intensity functions."^’* 


III. Forwardscattering Corrections 

A. Exact Mie Theory Formuiation 

^ mentioned earlier, because of the detector’s finite 
field of view, some forwardscattered light invariably 
wters into the detector’s view cone along with the direct 
beam. Bouguer’s law cannot be, strictly speaking, used 
to obtain the true optical depth r from transmission 
measurements. What one really obtains is the apparent 
opti^ depth t' related to the apparent efficiency fac- ' 
tors Q', which are distinguished from the true quantities 
Q; however, Q,b. ia clearly not affected by forward- 
scattering. 

■ Qmat + Qib* (apparent), ( 7 ) 


oerN 

orrtCToa 



<?Mt - QUt - Qmm - f*(ii + it) tiiO'di’ (8) 

X* %/o 

Note that for the small values of 9 considered in this 
work {9 < IQO), at and we may thus neglect the ef- 
fects of polarization. Then a correction factor due to 
forwardscattering can be defined as 

^ / Qnt ■ l - £, (9) 

where E, the error factor, is 

(10) ■ 

B. Approximate Rayleigh Diffraction Formulation 

Gumprecht and Sleipcevich^-2 have shown for small 
values of 9 and large values of x that as a good approxi- 
mation to the exact Mie theory the value of R can be 
computed from Rayleigh’s diffraction theory formula- 
tion9 for large particles. Rayleigh has shown by dif- 
fraction thTOry that the fraction of the diffracted light 
falling outside a cone of half-angle 9 in the forward di- 
rection is given by 

fit) ■ (xS) + (11) 

so that the forwardscattering correction and error fac- 
tors are defined as 

R • R(x«) - i [1 + Jj(x 4 ) + (12) 

^ ■ J - Jt(x«)), (13) 

where Jq and J i are Bessel functions of the first kind 
and of orders zero and one, respectively, ^ is in radians 
and 18 small enough so that 9 a* sinfl. We assume Q « 
2.0 for large particles. 

The Rayleigh diffraction approximation has been 
investigated by a number of authors, including Hod- 
Unson and Greenleaves,^ who also examined the error 
due to forwardscattered light using a formula similar 
to (13), In addition, they examined the geometrical 
optics effects of external reflection and transmission 
plus refraction, but their results indicate that these 
contributions are negligible for 9 i 10*. 

IV. Corrected Transmission Law 


1 r*xncLi 


sounce 




orrtcTOd 


Fig. 1. Schematic illustrations of (a) an open detector transmisso- 
meter system, (b) a lens-pinhole detector transmissometer system, 
and (c) the projection of an off-axis particle on a section of the open- 
detector system. 


By inserting the correction factor R into the trans- 
mission equation, one can account for both the direct 
and forwardscattered radiation. Thus, 

/(X) ■ /o<X) exp (-t'(X)), ( 14) 

where 

' X (15) 

Since the correction factor depends also on the geometry 
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of the transmission experiment, examples of two ex* 
perimental designs are discussed in the next section. 

V. Experimental Design Considerations 

An extinction experiment essentially involves a 
source of radiation (S), a transmitter lens to transmit 
a collimated beam of known intensity Iq(\) and radius 
Ri, through the aerosol medium of interest, and a de- 
tector system to measure the intensity of the direct 
beam after it has traversed a thickness L inside the 
medium. The correction factor R depends on the de- 
sign of the measurement system, and, therefore, the 
design must be given careful considerations. In light 
of the discussions in the previous sections, two simple 
experimental systems are compared for their perfor- 
mance. Both use a collimated beam of radiation; they 
only differ in the design of their detector systems. One 
uses an open photodetector system [Fig. l(a)J and the 
other, a lens-pinhole system [Fig. 1(b)]. The former 
will be referred to as simply the open detector sys- 
tem. 


A. Open Detector System 

The experimental system, shown in Fig. 1(a), consists 
of a radiation source, 5, situated at the focal length of 
a transmitter lens Li, which transmits a parallel beam 
through an aerosol medium of thickness L, which is 
measured by an open photodetector of radius The 
detector subtends different angles 6 for scatterers lo- 
cated at different locations I inside the beam. Thus, the 
results differ for particles on and off the beam axis. 
These two cases are discussed as follows. 

i. On- Axis Particles 

For particles on the beeun axis, angle 6 is given by 

tanS - R^iL - /). (16) 

In this case, since 8 « 8(1), Eq. (15) can be written as 

r' • Lwr^Np(r)Q^{x,m)R(xM), (17) 

where the path averaged correction factor R for a ho- 
mogeneous aerosol medium is given by 

/?(xX)-7 r^/eix,tf(/)idi 

L «/0 

X w/l 

R{x,S) coMiASdS, (18) 

to “ tan”' (Rj/L). 

In the exact formulation, the integral in Eq. (18) is 
solved by numerical quadratures. However, if one 
wishes to get some idea of how R behaves as a function 
of jc and L, it might be worthwhile to obtain an analyt- 
ical solution of Eq. (18) under certain reasonable ap- 
proximations. If one makes the approximations that 
tand ae 8 and L » R2 in Eq. (16), Eq. (18) reduces to 

R(xX)-^ f* R(x,t}t-^dt (19) 

L Jr%/l 

Using Eq. (12) and letting y = R^x/L, Eq. (19) yields 


/* ■ /Uy) - j ^ 

+ |y* (J|(y) ■♦•./f(y)l (20) 
3 3 

When y -^Q,R reduces to 

/?o - 1 - 4y/3» y V8 + 0(y«), 

and when y -• ®, i? reduces to 

/I- j + (2/(3r)| cM2y. 


In some experimental situations, there may be a gap 
between the scattering medium and the detector. In 
this case, the upper limit of integration in Elq. (18) will 
be less than r/2. Clearly, the integral can again be done 
numerically, but we may also use Eq. (20) if we subtract 
off the con^bution to the integral due to the empty gap. 
For example, consider the scattering medium of length 
L, and assume that the detector is located at a distance 
fj beyond the edge of the medium. Then it may easily 
be seen that 


R(l) m R(l + 

L L 


( 21 ) 


2. Off-Axis Particles 

For particles off the beam axis, the situation is clearly 
more complicated, as the cone subtended by the de- 
tector is no longer coaxial with the incident radiation. 
Thus, one can no longer assume the 0 integration im- 
plicit in (9), 0 being the azimuth. This additional 
complexity is beyond the scope of this work. 

However, it is not too difficult to get some idea of the 
errors involved in ignoring this complication. Consider 
a particle at a location I along the beam and a radial 
distance 6 normal to the beam axis, as illustrated in the 
projection shown in Fig. 1(c). Here 5 < ^ 2 - Then one 
may determine both the average scattering angle 8 and 
the average projected distance 3, where 


tane ■ d/lL — 1) 

R| ■ 4* + d* + 26d COM, 

5 ■ - r 4(e) de, 

T Jo 

3 ■ - r d(e)de. 

r Jo 

Although 8 cannot be determined analytically, 3 can be 
expressed in terms of the complete elliptic integral of 
the second kind^®; 


r {(I - p* sin^e)'^ - p coMlde 

r e/0 


2Rt 

■ — £(p). 

r 


( 22 ) 


where p ■ 8/R2- 8 was evaluated numerically for a series 
of values of p and (L — D/R2, and it was found that, to 
a very good approximation in most cases, tand * 3/(L 
— 1). This approximation was found to be invalid only 
for the cases when 5 was large, and (L - l)/R2 was smali 
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Thus, under most circumstances, the following equation 
is valid, namely. 

ub9 - - 1)\. (23) 

When p • 0, £ ■ t/ 2, Eq. (23) reduces to Eq. (16). The 
tables of E may be used to infer typical magnitudes of 
errors in d, if Eq. (16) is used. Thus, when p “ 0.5, ? « 
0.934£^ then the error in tand is ^ 6.6% high; when p * 
1/V2, a * 0.860A2* tbe error is '>'14% high; and for p as 
high as 0.9, ? « 0.745i?2, so that the error is about 25% 
high. 

If particles are uniformly distributed out to some 
maximum value Pm»x of p, the average value of 1 for 
these particles is 

^ " Jo j updp 

- ^ K(P^] . (24) 

where K is the complete elliptic integral of the first 
Itind.^® 

One may now determine the average values of errors 
in ? and hence 9, involv^ in ignoring off-axis particles. 
Thus when p^ax * 0.5, a * 0.968£2> resulting in a 3.2% 
error, •when Pn»** * l/v'2,^ ■ 0.914/? 2 . an 8.6% error; 
even when Pn»*x * 1. then d ■ 0.849/?2, resulting in an 
error of only 15.1%. 


B. Lens-Pinhole System 

The experimental system illustrated in Fig. 1(b) 
consists of a radiation source, S, situated at the focal 
point of a lens L\, which produces a collimated beam of 
light. At the end of the beam path is a second lens L 2 
of focal length /, which focuses the light through an 
aperture of radius a in the focal plane and onto a de- 
tector. For such a system, it is easy to show that the 
half-angle 6 is constant and is given by 


VI. Computational Considerations 

A parametric study of the correction factor R and the 
corresponding percentage error £ was carried out as 
functions of different combinations of x, 9, m', and m', 
whose values occur within the ranges 0.1 ^ z ^ 100, 0* 
5 fl i 10®, 1.33 5 m' S 1.65, and 0.0 i m' i 0.1. The 
somewhat high values of 9 and m', such as » 10* and 
m' ■ 0.1, were intentionally selected in order to un- 
derstand the behavior of R and E aa 9 and m' approach 
these large values to gauge how valid are the approxi- 
mations made in the theory in that region. 

In the case of computations for R (and £) based on 
Mie theory, instead of employing a numerical quadra- 
ture scheme for evaluating the integral in Eq. (10), 
namely. 


/(x,m,S) • r (ii it) •int'df' ■ “ T* (»i iiidit', (27) 

Jo Jo 

where p « cosd and p' « cos$', use was made of the fol- 
lowing explicit closed form solution given by Wiscombe 
and Chylek [Eq. (4) in Ref. 11] 


- (1 - M*) U1/Ii(m)| + M|il2(M)| - 2 R»(«i(p)il5(M))l 


where 


+ i (2> + l)(K|»+ |6k|*)/f*(M) 
*•1 

• R>(dtd(* -f 6>6|*) 
nmt (•! /({ -f i) — k{k + 1) 

X [P|(M)rk(p) - Pk(M)«l(M)l(l - M*), 


Hk 

dk 

Ri(m) 


- (2* + 1) J^"pfU')dM'. 

- (2* + - (2* + Dbk, 


• dkTi(ll) 
*-!*(*+ 1) 




• Skruiit) I 

*-i *(*+1). / 


(28) 


(29) 


and where Pi, is the Legendre polynomial, /?2 represents 
the complex conjugate of R 2 , and a*, bk, and vi, have 
their usual Mie theory meanings."*^ Then 


Uni) <■ a/f (25a) 

provided the following condition holds, namely, that 
RtiRi, + aUf. (2Sb) 

where /?2 is the radius of lens L 2 . Since 9 is constant for 
all particles anywhere inside the beam, for this system, 
referred to as the lens-pinhole system, Eq. (17) can be 
written as 

r' • Lrr»NpQ„,(x.m)R(x.S). (26a) 

since 

R - Rix.S), (26b) 

as given by Eqs. (9) or (12). 

For particles whose distance from the beam axis 5 
and distance from the detector d jje such that 5 + d- a/f 
> R 2 , the collecting lens L 2 may be considered a flat 
detector, and the results obtained earlier employed. 
Thus for circumstances in which Eq. (25b) does not 
hold, it is necessary to employ a combination of re- 
sults. 


f:-/(i,m.»)/**<?«. (30) 

Computations for R and £ by the approximate method 
were made by using Eqs. (12) and (13). Computations 
for the path-averaged correction factor R were made for 
various values of y by using Elqs. (12) and (20) for the 
lens-pinhole and open detector systems, respectively. 

In Eq. (18), /? is a function of the two variables z and 
^Of which makes their plotting too complicated. In Eq. 
(20), however, /? is a function of y only, which may be 
far more easily visualized. Equation (20) is, of course, 
only an approximation to Eq. (18), and it is important 
to determine the errors involved in using this expres- 
sion. For this reason, Eq. (18) has been evaluated nu- 
merically for a series of values of z and 9q, and the re- 
sults have been compared with those of Eq. (20). 

The largest discrepancies occurred for large 9q and for 
z * 1.0 to 2.0. For smaller values of z, both results were 
quite close to 1.0 for 9q < 20*, although the errors in (1 
— /?)-were fairly large. Note that the use of the Ray- 
leigh expression in the evaluation of Eq. (18) probably 
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Fif. 2. Plots of the correction factor R vs Fig. 3. Plots of the percentage error (lOOc) Fig. 4. Plots of the percentage error (100£) 
size parameter x for refractive -index values: as a function of the size parameter z for the as a function of the size parameter x for re* 
1.33, 1.65, 1.50, 1.50-<(0.02); and 1.50- refractive-isKlez valuae: 1.33, 1.65, and fractive-indez values: 1.50, lJO-t(0.02), and 

»(0.1). 1.50. 1.50-«(0.1). 


introduces as much error as the approximation tand at 
6 for X values in this range. 

For ^0 * 20", the largest discrepancy was 10.5% when 
X ■ 2. For 6o * 10", the largest discrepancy was 4.6%, 
while for Bq « 5", the largest discrepancy was only 2.2%. 
For smaller $o values, the errors decreed proportion- 
ally. 

During this study, two other calculations were made. 
First, X sin^ was us^ instead of xB as the argument for 
the Bessel functions in the Rayleigh expression for R. 
This produced virtually no effect on the values of Eq. 
(18). Second, the integral in Eq. (19) was terminated 
at t/ 2, instead of «. This was found to produce con- 
siderable errors in almost all cases studied. 

VII. Discussions and Conclusions 

As mentioned earlier the parametric study of R was 
carried out as a function of many different combinations 
of X, B, m', and m", but for the ^e of clarity the results 
of only a few judiciously selected combinations are 
presented in the following subsections. 

A. Results for R and E as Functions of x and m 

Figures 2, 3, and 4 graphically depict the behavior of 
R and E as functions of x in the 0.1-100 range, for (1) 
three typical values of the view cone half-angle B, 
namely, 1", 4®, and 10"; (2) three typical values of m', 
namely, 1.33, 1.50, and 1.65; and (3) two complex values 
of m, namely, 1.50-t0.02 and 1.50-t0.1. In reference to 
the plots E denotes percentage error obtained by 
multiplying E by 100; elsewhere, it denotes simply the 
error factor E. 


For smaller x. Fig. 2 shows that for the three values 
of 8 and the five values of the reffactive index, the value 
of R remains close to unity so that many details of the 
plots cannot be discerned. On the other hand, for the 
same small values of x the semilog plots of E vs x (Figs. 
3 and 4) depict the behavior of E in much greater detail 
The small'particle limit behavior of E can be easily ex- 
pressed by approximate analytic expressions as follows. 
As X values approach the small Rayleigh particle limit, 
one can use the lowest order approximations for I'l, ij, 
®nd Qtn given in Sect. 3.9 of ^f. 12 to derive the fol- 
lowing approximate expressions for E for both the 
nonabsorbing and the absorbing particles, namely. 


E{m’ « 0) - 'li - y* CO90 - % co»^9 ai 3a*/8, (31«) 


£(m' 0) «• £(m' 


0 ) 


L/n' 9m' J 


9*x^ (m'« - !)» 
m' 24m' 


(m' <C 1). 


(31b) 


Table I presents some sample approximate results, 
calculated with Eq, (31a) and (31b), along with the exact 
results, for ^ » 1* and 10", x » 0.1, 0.5, and 1.0, and m 
» 1.5-i(0.0), 1.5-i(0.02), and 1.5-i(0.1). From these 
results one can see that the approximate formulas are 
good to within a few percent of the exact results for x ■ 
0.1 and are still reasonably good for x » 0.5. However, 
for X > 1.0, the behavior of E differs considerably from 
that predicted by Eq. (31), so much so that the m" de- 
pendence virtually disappears for x > 3.0. It should ^ 
noted, however, that the 6^ behavior, predicted by Eq. 
(32) for both absorbing and nonabsorbing particles, is 
Quite valid for 6 < 10® and x 5 5.0. 
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For large values of x (x ^ 10) the effects of the variojM 
diffraction maxima and minima, or corona,® can be 
clearly seen, particularly in Fig. 2. These effe^ ut 
more important than the effects of refractive index for 
particles in this range, although it should be noted that 
the presence of some absorption tends to damp the os- 
cillations. It should be pointed out that ^y attempt 
to repeat these graphs with a higher density of pomts 
in the 10 i X 100 cycle merely results in a fm more 
cluttered web. These oscillations disappear when the 
intensity is averaged over a suitable size distribu- 
tion.^® _ , , 

Figure 3 represents the behavior of E for noimb- 
sorbing particles, and Fig. 4 represents that for partcles 
with different values of the imaginary part of the re- 
fractive index. For nonabsorbing particles, the values 
of £ rapidly reduce to constant threshold values 
(.^3^/8) independent of x and m', as x reduces to less 

than about 1.0. . ^ r j n 

For absorbing particles, the behavior of £ and a re- 
main almost identical to that of nonabsorbing particlM 
for this range of tf and x > 1.5. However, for x ^ 1.5, £ 
decreases rapidly to values much lower ti^ the 
threshold value for nonabsorbing particles. For ve^ 
.mall absorbing particles the absorption dominates the 
scattering, £ varying inversely as m'. 

B. Results for R and B as Functions of 6 
Figures 5 and 6 illustrate the behavior of R and £, 
respectively, as functions of B in the range 0.25“ to 10 
for m - 1.50 and for eight different values of x from 0.1 
to 100. In Fig. 5, one sees that for x ^ 5.0, the valu^f 
R falls off quite rapidly from unity, as B increases, ^e 
plots of £ vs in Fig. 6 depict the percentage errors due 
to forwardscattering for small values of tf m greater 

detail In Fig. 6, one sees that £ increases as Omcreases 

for all values of x, until it reaches a maximi^ value of 
about 0.5. The plots for x - 0.1 could easily be repre^ 
sented by Eqs. (31). 


T«M*L Coni|»«rt«»n ^ 

lor fho tiTor Foctor 


Rrr or factor 

X TppFSr Eisct Aporoa. BiMrt 


m 

1.5-^(0.0) 

X 

0.1 

0.5 

1.0 

1.142(-4) 

1.142(-4) 

1.142(-4) 

1.148(-4) 

1.278(-4) 

1.738(-4) 

1.142(-2) 

1.142(-2) 

1.142(-2) 

1.138(-2) 

1.266(-2) 

1.716(-2) 

1.54(0.02) 

O.l 

0.5 

1.0 

6.610(-7) 

8.283(-5) 

6.610(-4) 

6.582(-7) 

5.012(-5) 

l.372(-4) 

6.610(-5) 

3.263(-3) 

6.610(-2) 

6.5l6(-5) 

4.958(-3) 

1.355(-2) 

1.54(0.1) 

0.1 

0.5 

1.0 

1.322(-7) 

1.653(-5) 

1.322(-4) 

1.375(-7) 

1.51K-5) 

0.76K-4) 

1.322(-5) 

1.653(-3) 

1.3221-2) 

1.38K-5) 

1.495(-3) 

0.75K-2) 


C. Error Contour Diagram 

It is obvious from the results shown in Figs. 2-6 that 
one can plot a series of percentage error contours in the 
X ~B plane to illustrate trade-offs that can be made 
between the values of x and B. Three contoim for 
percentage errors of 1%. 5%, and 10% are shown m Fig. 
7 for m » 1.5-i(0.0) and for the ranges 1 $ x £ 100 and 
0“ :S ^ 10“. T^e 1% contour shows an inflection pomt 
in its curvature at about d a 5". 



Fi*. 5. PloU of the comction factor R va half-angla • for m “ 1.50 
and eight values of the siza parameter x. 



Fig. 6. Plots of the percentage error ( 100£) vs half-angle for m ■ 
1.50 and eight values of the size parameter x. 
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Fi|. 7. Error contour* inx —$ pUn* for m ■ 1.50 and pcrcenUg* 
•rrorc 1%, 5%, and 10%. 




Fif . 8. PloU of A va half-angle 8 for m » 1.33 and five value* of x. 
(Symbol* without croase* represent Mie theory results; those with 
crusses, approximate theory results.) 


D. Comparison Between Mie Theory and Diffraction 
Formuiations 

The approximate expressions for R and E in Eq. (13) 
have no m -dependence. Figures ft-12 show plots of R 
vs as computed using both Eqs. (11a) and (13a), for 
a series of vdues of * and a series of refractive indices. 
Symbols without a cross represent exact Mie results 
[Exj. (11a)), and symbols with a cross represent Rayleigh 
diffraction results [Eq. (13a)]. 

Figures 8, 9, and 10 show that for nonabsorbing peur- 
ticles with three different refractive indices (1.33, 1.65, 
and 1.5), the Rayleigh diffraction results are in fairly 
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Fif. 9. Plott of A as a function of 6 for m * 1.65 a.nd fiv# valuM of 
X. (Symbota without croaaas rapresont Mi# thaory raauJta: th«M with 
croas#*, approzimat# theory results.) 


Fif. 10. Plot* of A as a function of 6 for m ■ 1.50 and five values of 
X. I Symbol* without croaaes represent Mie theory results; those with 
crosses, approximate theory resulto.) 


good agreement (within about ± 5%) with the exact Mie 
results for a wide range of x and 9 values. However, for 
small values of x, while both expressions yield results 
close to unity for R, they may yield significantly dif- 
ferent values of E. The agreement appears to be best 
for 9 < 1.5". 
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Fi«. 11. Plou of A at c function of • for m • l.50-i(0.02) and fiva 
valuM of X. (Symbols without croasaa rtpretant Mia thaory raaulta; 
thoaa with croasaa. approzimata thaory rasulu.) 


For absorbing aerosol particles. Figs. 11 and 12 show 
that agreement improves rather dramatically for 8 larger 
than 1.5* as m' increases its value from 0.00 to 0.1. 
Even when the absorption is small, the agreement be- 
tween results of the two methods significantly improves, 
especially for large values of x. This is an extremely 
useful result when absorbing aerosol particles are in- 
volved, as the use of the Rayleigh diffraction formula 
[Eq. (13a)] can yield accurate results, while at the same 
time considerably reducing computation costs. In 
addition, the approximation might be amenable to an- 
alytical integration in some cases, as shown earlier, 
which could be a valuable asset in the optimum de- 
signing of an extinction experiment. 

E. Optimum Transmissometor Design 

Equation (14) shows how one includes the for- 
wardscattering correction factor to obtain the measured 
(or apparent) optical depth r'. Of the various optical 
systems that have been used in light-transmission ex- 
periments, two of which are described in eariier sections, 
the lens-pinhole optical system is preferred because the 
value of is constant and independent of the location 
of the illuminated particle inside the beam path, as long 
as Eq. (25b) is satisfied. In addition, the lens-pinhole 
system practically excludes all stray light from the de- 
tector. 

In contrast to the lens-pinhole system, the open de- 
tector system {Fig, l(a)j presents obvious difHculties, 
such as (1) ^ is not constant but depends on the location 
of the scatterer in the beam path, so that computations 
for H are somewhat cumbersome: and (2) the value of 
R for particles near the outer edges of the beam differs 


Fif. 12. Plou of A as ■ function of f for m “ 1.50-i(0.1) and fiva 
valuos of X. (Symbols without croasaa rapraaant Mis thaory rasulU; 
thoaa arith croasaa. approzimsu thaory raaulta.) 


from that for particles at the center of the beam at any 
given distance from the detector, due to the fact that the 
axis of the cone subtended by the detector with the 
former particle is not coincident with the direction of 
propagation of the incident light. 

The values of the path averaged correction factor R 
as given by Eqs. (20) and (26b) for the open and the 
lens-pinhole detector systems, respectively, are plotted 
in Fig. 6 as a function of y (» x8), where R 2 /L and 
tan'^a// are the respective vsdues of 8 for the two sys- 
tems. The values of R for the open detector are smaller 
than those for the lens-pinhole system for the 0.01 < y 

< 100 range. This difference in the two sets of values 
of R tends to disappear for extreme values of y, viz., y 

< 0.02 and y > 15.00. For y — «, /? — 0.50. These 
results are useful for designing transmission experi- 
ments. Maximum accuracy is achieved if one designs 
the experiment such that R is nearly constant, which 
occurs at the extreme values of y . In these y ranges, the 
lens-pinhole detector ‘s more accurate than the open 
det^tor (see Fig. 6). However, in the intermediate y 
region, the uncertainty in R due to variations in y for the 
lens-pinhole detector can be larger than that for the 
open detector due to the steeper gradient. 

From the foregoing discussion, it is apparent that the 
transmission law as ordinarily expressed by Eq. (1) 
caimot be used except in the case of small particles, for 
which it is only a good approximation. The effect of not 
only the particle’s size and composition but also the 
geometry of the optical system on the apparent volume 
scattering and extinction coefficients must be consid- 
ered. 

The results discussed in this paper pertain to the case 
of a collimated beam of radiation traversing a homo- 
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Fig. 13. PloU of th« path-avtrag od comction factor /} as a function 
of y, whara y • xa/f and xR^/L for tha lana-pinhola datactor and 
opan-datactor syatama, raapactivaly. 

geneous monodisperse aero«ol medium under the con* 
^tions in which the single-scattering approximation is 
valid. The case of the homogeneous polydisperse 
aerosols is presented in a separate paper (Part 2) in this 
journal Recent experimental investigations^^ of the 
attenuation of argon laser beam in fogs (with different 
size distributions of droplets) as a function of the re- 
ceiver field of view (between 0.5” and 4.0*) confirms our 
theoretical model for R as a function of 9 . 


This work was performed under NASA grant 1252 
and contract NASl*15198. 


References 

1. R. 0. Gumprvcht and C. M. Sliapctvich, J. Phya. Cham. 57, 90 
(1953). 

2. R. 0. (Sumpracht and C. M. Sliapcavich, J. Phya. (Tham 57, 95 
(1953). 

3. A. Oaapak and 0. H. Vaughan, AppL Opt. 17 , 374 (1978). 

4. J. R. Hodkinaon and L Graanlaavaa, J. Opt. Soc. Am. 53, 577 
(1963). 

5. G. E. Shaw, Pagtoph (BirkhSuaarVarlag, Baaai, 1976), p. 114. 

6. R. W, AuadnandT.J. Patzoid,‘‘ConaidanaonaintfaaDaaignand 
Evaluation of O»anographk Tranaitiimomatara,** in Light in the 
Sxa. J. E. Tylar, Ed. (Acadamic Praaa Naw York, 1977). 

7. R C. van da Hubt, Scatttring of Light hy SrtuUl Partieln (Wilay, 
Naw York, 1967). 

8. 0. Dairmandjian, EUetromagnttie Scatttring on Sphtrical Po~ 
lyditptrtiont (Elaaviar, Naw York, 1969). 

9. Lord Raylaigh, Philoa. Mag. Sar. 5 12, 81 (1881). 

10. L S. Gradahtayn and L W. Ryzhik, Tablet of Integralt, Sariet and 
Produett (Andamic, Naw York, 1965). 

11. W. J. Wiacomba and P. Chylak, J. Opt Soc Am., 67, 572 
(1977). 

12. M. Karkar, The Scattering of Light and Other Electromagnetic 
Radiation (Acadamic, Naw York, 1969). 

13. A. Daapak and M. A. Box. “Forwardacattaring Corractiona for 
Optical Extinction Maaauremanta in Aaroaol Madia. 2: Poly* 
diaparaiona” 17, 0000(1 Oct., 1978). 

14. A. Carawall, York U., Toronto; private communication. 


34 




APPENDIX 2 


FORWARDSCATTERING CORRECTIONS FOR OPTICAL EXTINCTION 
MEASUREMENTS IN AEROSOL MEDIA; 2: POLYDISPERSIONS 


Reprintad from APPLIED OPTICS. VoL 17, pace 3169, October 1, 1978 
Copyrifht 1978 by the Optical Society of America and reprintM by permisaion of the copyright owner. 


Forwardscattering corrections for opticai extinction 
measurements in aerosol media. 2: Polydispersions 


Adarsh Deepak and Michael A. Box 


ORIGINAL PAGE IS 
OF POOR QUALITY 


This paper, second of two parts, presents a parametric study of the forwardscattering corrections for experi- 
mentally measured optical extinction coefficients in polydispetee particulau media, since some forward 
scattered light invariably enters, along with the direct beam, into the finite aperture of the detector. For- 
wardscattering corrections are computed by two methods: (1) using the exact Mie theory, and (2) the ap- 
proximate Rayleigh diffraction formula for spherical particles. A parametric study of the dependence of 
the corrections on mode radii, real and imaginary parts of the complex refractive ind*x, and half-angle of the 
detector’s view cone has been carried out for three different size distribution functions of the modified 
Gamma type. In addition, a study has been carried out to investigate the range of these parameters in which 
the approximate formulation is valid. The agreement is especially good for small-view cone angles and large 
particles, which improves significantly for slightly absorbing aerosol particles. Also discussed is the depen- 
dence of these corrections on the experimental design of the transmissometer systems. 


I. Introduction 

In Part of two papei^ti, results were presented of 
a parametric study of the forwardscattering correction 
factor R and the complementary error factor E for 
monodispersions. In t^ paper the results of a similar 
study carried out for spherical polydispersions of size 
distribution n(r) (cm~® r being the radius in ^m, 

will be presented. For the sake of clarity, only the re- 
sults obtained with the use of simple unimc^al size 
distributions of the modified Gamma type, such as 
Dermendjian models^ Haze M, Haze H, and Cloud C3 
(referred to as Haze C in this paper), are presented here. 
Results for other real size distributions can easily be 
obtained in a similar manner. The behavior of both R 
and E, averaged over each of the three size distributions, 
will be discussed here as functions of each of the fol- 
lowing parameters: the mode radius r„; the piolydis- 
persity or the spread of the size distribution; and the real 
(mO and imaginary (m') parts of the complex refractive 
index m ■ m' - im*. The computations have been 
carried out with both the exact Mie theory solutions in 
explicit closed form and the Rayleigh diffraction theory 
approximation, as explained in Part 1. The results of 
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such a study are extremely useful in obtaining the op- 
timum experimental design parameters for the mea- 
surement of extinction coefficients in particulate 
media.^ 

II. Transmission Law (Bouguer’s Law) 

The transmission law for an electromagnetic plane 
wave passing through a homogeneous poiydisperse 
aerosol medium (Fig. 1, Part 1) is given by Bouguer’s 
law, defined by Eq. (1) in Part 1, namely, 

/fX) ■ /o(^) sxp(-T<X)l, . (1) 

where 

HX) - d«(X)L. (2) 

dm(X) - -(■ d.ta(X). (3) 

As in Part 1, for the sake of clarity, the molecular con- 
tributions to be ignored here. Then for poly- 

disperse aerosols of size distribution n(r), cm~^ Mm~‘, 
the coefficients are defined by 

dy(X) ■ r »T*Q;(x,m)n(r)dr,; • ext, scat, and ab». (4) 

All the quantities in Eq. (4) are the same as defined in 
Part 1, except that the factor iVp(r) for monodisperse 
particles has been replaced here by the operator 



where ri and rj are lower and upper limits of radii. 
Aerosol size distributions are discussed in a later sec- 
tion. 
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/(X) - /o(X) •xp(-r'(X)), 


( 11 ) 


III. Forwardscattaring Corractlons 

A. Exact Mie Theory Formulation 

Because of the fact that forwardscattered light in> 
variably enters the detector view cone, Bouguer’s law, 
as defined in Eq. (1), cannot be used to obtain the true 
optical depth r in a transmission experiment. But in* 
stead, one obtains the apparent optical depth r' related 
to the apparent volume extinction scattering 

coefficients, which are distinguished from the true 
quantities and d«cMi respectively, here, dab* i* clearly 
not affected by the forwardscattering. T^us, 

” (fiM (apparant). (5) 

The Mie efficiency factor is based on the total 
amount of light lost from the beam by scattering by a 
particle in all directions; and Q\c»t> on the amcunt of 
light lost by scattering in all directions except within a 
cone of half-angle 9 in the forward direction, so that 
their difference is given by 

Qt*i ” Qax “ Qtctt ” Qmm " “I f* Ui ^ (i) ainO'dS 
X* Jo 


where I(x,m,d) is defined in Eqs. (28) and (29) in Part 
1. From Eqs. (3) and (5), one obtains 

~ dut ■ <^ic« ■“ ^km 

■ » '•»fi(r)((}^ - Q^dr 

J'l 

■ *• r^(r)QEdr 

Jr\ 

C'* n(r)l{x,m,$)dr, (7) 

** Jrt 

where E (and R) are defined in Eqs. (10) and (9) of Part 
1. Then a correction factor R and an error factor E, 
averaged over the particle size distribution between the 


limits and r 2 , may be defined by 


R • • l - F. 

(8) 

where 


R ■ C'' n(r)Hx,m,B)dr. 

(9) 

-In 


Note that R ■ ff[n(r),m,tfl. 



B. Approximate Rayleigh Diffraction Formulation 

From Eq. (12) in Part 1 and Eqs. (6) and (8), the 
Rayleigh diffraction approximation to the forwards- 
cattering correction for a polydispersion is given by 

ff - ~ f * r*n(DQ(x.m)‘<,(l + JUxB> + J^xt)\dr, (10) 

where Jq and </i are Bessel functions of the first kind 
and of orders zero and one, respectively. 

IV. Corrected Transmission Law 

By inserting the correction factor R into the trans- 
mission equation, one can account for both the direct 
and forwardscattered radiation. Thus 


where 

t'(X) - 

From the discussion in Part 1, it is clear that for parallel 
beam transmission systems with an open detector, the 
half-cone angle 9 * 9(1), so that 

r'(X) • tL r''dfTHJ(x./n)fi(r)/l(xX), (13) 

Jn 

where 

A(xX )-7 f‘dW(x,»(01. (14) 

L vO 

And for a lens-pinhole system 9 >■ constant, so that R 
•R(x,9), independent of /. For the discussion of the 
experimental design considerations of the two detector 
systems, see Part 1. 

V. Aerosol Size Distributions 

Several analytic representations of aerosol size dis- 
tributions appear in the literature. In this paper, the 
behavior of R and E is investigated as a function of the 
mode radius r„, the spread of the size distribution and 
the complex ref^ctive index m' im"), the upper 
and lower limits of radii being 10"^ and 20 4m. For the 
sake of simplicity, three of the Oeirmendjian models,^ 
namely. Haze M, and Cloud C3 (referred to here as Haze 
C), were selected for representing different polydis- 
persities of aerosol size distribution n(r). Since is 
varied between the radii limits and 10 iim, the three 
models are used here in a more general way than was 
their original intent. It is in that sense that model 
Cloud C3 is referred to as Haze C. The limits of inte- 
gration over r are 10“* and 20 4m. 

The expressions for the n(r) models and their corre- 
sponding mode radii r„ are given as follows: 


(a) Haze M: 

n(r) ■ r *xp|-(6r)‘^*], r„ • 4/b, 

(15) 

(b) Haze H: 

n(r) ■ r* exp(-6r), - 2/6, 

(16) 

(c) Haze C: 

fi(r) ■ ’■•exp(-(6r)*l, ■ (8/3)‘^®6"‘. 

(17) 


Note that any normalization or scale factor in the size 
distribution will cancel when E and R are evaluated. 
Thus, although E and R do depend on the shape of the 
size distribution, they do not depend on the total 
number of particles. 

The difference in the shape of these three models is 
illustrated in Fig. 1, where each has r„ • 1.0 4m and 
n(r„) ■ 1.0 cm“ For want of a better terminol- 

ogy, the spread of a size distribution will also be referred 
to as polydispersity of the size distribution in this paper. 
For example. Haze H will be referred to as more poly- 
disperse than Haze C, and in the same vein it will be 
stated that Haze M has a higher polydispersity than 
Haze H. It should be noted here that there exist other 
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Fig. 1. Plot* of the three Deirmendjien size-distribution models. 
Hazes M, H, and C, for a mode radius of 1 Mm. 


tive region. For this reason, we have also considered a 
power law distribution, given by 

n(r) ■ r“», 10“* Mm S r 2 15 Mm 

and allowed r to vary between 2 and 4. Note that one 
cannot talk about either mode radius or polydispersity 
for a power law haze, only slope i>. 

VI. Computational Considerations 

A parametric study of the correction factor 7? and the 
corresponding error factor E was carried out as func* 
tions of different combinations of r„, B, m' and m", 
whose values occur within the ranges 0.1 < < 10 Min, 

0^ <B< lO", 1.33 <m'< 1.65, and 0.0 < m' < 0.1. 

The computations of (and E) in Eq. (8) [and Ek). 
(9)] are made by using the closed form relations given 
in Part 1 for I(x,m,B). The computations of E for the 
approximate method are made by using Eq. (10). 

VII. Discussions and Conclusions 

The parametric study of E was carried out as a 
function of many different combinations of r„, B, m', 
and m', but for the sake of clarity the results of only a 
few judiciously selected combinations are presented in 
the following subsections. As the Deirmendjian models 
are so different from the power law model, we shall treat 
each separately, starting with the Deirmendjian mod* 
els. 


terms in the theory of distribution functions to express 
the same quantity. 

Perhaps the major disadvantage of the Deirmendjian 
models is that their falloff behavior for large radii is too 
sharp; many experimentally measured distributions 
show a power law behavior, at least in the optically ac- 


A. Results for R and E as Functions of r„ and m 

Figures 2-7 illustrate the behavior of E and E, re- 
spectively, as functions of r„ in the 0.1-10.0-Min range, 
for three different half -cone angles (1*, 4®, and 10*) and 
five different refractive indices [1.33, 1.65, 1.55, 1.55- 
i(.05), 1.55-{(0.1)], for each of the three models: (a) 



Fig. 2. Plots of ths correction factor vs 
mode radius for three values of 9 and five re- 
fractive indices: Haze .Vf . 



Fig. 3. Plots of the correction factor ff vs 
mode radius for three values of 9 and five re- 
fractive indices: Haze H. 



■MODE SAOUS „Mn 


Fig. 4. Plots of the correction factor Tf vs 
mode radius for three values of 9 and five re- 
fractive indices: Haze C. 
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Fig. 5. Plot* of percentage error (lOOE) V8 
mode radius for three values of 9 and five re- 
fractive indices: Haze M. 


Fig. 6. Piou of percentage error (lOOf) vs Fig. 7. PloU of percentage error (lOOff) vs 
mode radius for three values of t and five re- mode radius for three values of 9 and five re- 
fractive indices: Haze H. Tractive indices: Haze C. 


Haze Af, (b) Haze H, and (c) Haze C, respectively. The 
computations were made for X ■ 0.56 nm. 

For Haze M (Fig. 2), the values of rapidly decrease 
as r„ increases from 0.1 ^m to about 1 tun and then tend 
to level off to a nearly constant value. The percentage 
error (Fig. 5) for Haze M shows that for small r„ the 
error E increases fairly rapidly and then quickly levels 
off to a constant value, which is different for each set of 

and m values. The smaller the d the higher is the 
value of Tm beyond which the leveling of the values for 
E takes place. 

For Haze H and Haze C, the plots of make an in- 
verted integral sign ( 2 ) within the r„ range of 0 . 1 - 10.0 
fim. their steepness increasing with the decrease in po- 
lycUspersity of the size distribution. From Figs. 2-7, we 
see that for ^ ■ 1 ®, there is virtually no m -dependence 
for any of the hazes. For this reason, very few symbols 
have been drawn on these lines. For S ■ 4®, we see that 
a gTWflll m -dependence has started to appear, being most 
clearly visible in Figs. 5-7, due to the logarithmic scale 
for E. However, for » 10 ® , a clew m -dependence cm 
be distinguished, especially in Figs. 2-4. We see, in 
general, that is lowest for m 1.33, and that there is 
little difference between the curves for m * 1.55 and 
1.65. However, the presence of a small amount of ab- 
sorption immediately raises E, although the actual value 
of this absorption appears to have little effert. 

Effect of polydispersity: Comparing the R (and E) 
vs r„ plots for the three models, one can reach the fol- 
lowing conclusions: 

It can be seen from the plots presented here that, in 
general, for given values of B and rm, the higher the po- 
lydispersity the higher the error E. This is due to the 
predominant forwardscattering of the large particle 


component, implied by the increased polydispersity. 
We should note, however, that the saturation values of 
E and R tboth 0 . 5 ) are not affected by polydispersity, 
but that these values are reached sooner for a more 
polydisperse haze than for a less polydisperse haze. A 
saturation value o(E ^ 0.5 implies a Q„t value of 1 . 0 . 
An explanation for this is as follows. F or the same value 
of B, the amount of scattered radiation collected at the 
detector increases as the large particle component in- 
creases, which, in effect, reduces Q' from a value of 
about 5 L 0 (assumed for large particles) to a minimum 
value of 1 . 0 . The saturation value for E tends to be 
about 50% for B « 10® and 45% for * 4® for Haze Af , a 
difference of about 5%. This difference, however, tends 
to decrease as aerosol size distributions become less 
polydisperse, implying a smaller number of large par- 
ticles, so that the satviration value of E approaches 
50%. 

B. Results for R and £ as a Function of B 

Figures 8-10 show the behavior of E as a function of 
B for six values of mode radius (0.1, 0.15, 0.50, 1.0, 2.5, 
and 10.0 ^m), and m ■ 1.56 for each of the three size 
distribution models, respectively. Again, the effects 
of polydispersity are quite apparent in the increasing 
spread of the E vs plots for lower r„ values (below 1.0 
4 ) as the polydispersity decreases. 

Figures 11-13 illustrate the behavior of ff as a func- 
tion of B for five values of r„ (0.1 4 m. 0.4 4 m, 1.0 4 m, 4.0 
4 m, and 1.0 4 m). Symbols without crosses represent 
the Mie theory results, and those with cross represent 
the approximate results. The discussion of their 
comparison will be presented in the next section. Only 
the Mie theory results will be discussed here. The 
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Fif. 8. Plots of percentage error (100?) ve 
half-angle 9 for six mode radii (m “ 1.55): 
HazeAf. 



«Moa. oca 


Fig. 9. Plots of percentage error (100?) vs 
half-angle 9 for six mode radii (m » 1.55): 
Haze H. 



Fig. 10. Plots of percentage error (100?) vs 
half-angle 9 for six mode radii (m - 1.55): 
Haze C. 





Fig. 1 1. Plots of correction factor ? vs 0 for 
m ■ 1.55, \ > 0.55 um. and five mode radii: 
Haze S4. I Symbols without crosses represent 
Mie results: symbols with crosses, approxi- 
mate results.) 


Fig. 12. Plots of correction factor ? vs for 
m * 1.55, X * 0.55 am. and five mode radii: 
Haze H. < Symbols without crosses represent 
Mie results: symbols with crosses, approxi- 
mate results.) 


Fig. 13. Plots of correction factor ? vs 4 for 
m > 1.55, X « 0.55 am, and five mode radii: 
Haze C. (Symbols wivhout crosses represent 
Mie results: symbols with crosses, approxi- 
mate results.) 
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Fig. 14. Plot* of correction factor Tfvt${or 
n ■ 1.55-1 (0.05), X • 0.55 ma, and five mode 
radiu Haze M. (Symbols without crosses 
represent Mie results; symbols with crosses, 
approzimau results.) 


Fig. 15. Plots of correction factor 7? vs for 
A ■ 1.55-t(0.05), X ■ 0.55 iim, and five mode 
radii: Haze H. (Symbols without crosses 
represent Mie results; symbols with cross, 
approximate results.) 


Fig. 16. Plots of correction factor 7? vs for 
m • 1.55-t(0.05), X ■ 0.55 Mm, and five mode 
radii: Haze C. (Symbols withoit crosses 
represent Mie results; symbols with cross, 
approzimau results.) 


choice of different r„ values for the cases shown in Figs. 
8-13 was merely for the sake of the clarity of the 
graphical display. The plots in Figs. 11-13 also show 
that the correction factor ? falls rapidly with increasing 
9 for large r,„ values and then levels off quickly to a value 
close to 0.5. For large particles, the plots are practically 
identical for all the three models. However, for small 
values of there is a considerable spread in the values 
of??. 

C. Comparison Between the Mie Theory and 
Diffraction Formulation 

Consider Fi^. 11-16 representing the ?7 vs ^ results 
for nonabsorbing (m « 1.55) and absorbing (m * 
1.55-10.05) aerosols, respectively. Symbols with crosses 
represent the results due to Rayleigh formulation and 
those without crosses the results due to Mie theory. 
Figures 11 and 14, 12 and 15, and 13 and 16 represent 
the models Haze M, Haze H, and Haze Cj respective- 
ly- , 

Figures 11-13 show that, for nonabsorbing aerosols, 
the agreement between the Mie and the Rayleigh ap- 
proximation is not as good for small r„, as it is for larger 
r„ for all three models. This is to be expected since the 
Rayleigh diffraction formula Eq. (13) is only valid for 
large particles. On the other hand, the agreement for 
both large and small To, appears to be better for Haze 
M (more polydisperse aerosol) than for Haze C, indi- 
cating that for more polydisperse size dispersions the 
agreement between the Rayleigh formulation and the 
exact formulation tends to improve. Comparison of 
these results with those for monodispersions in Part 1 
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also substantiates the trend. The remarks for the cas e 
of nonabsorbing aerosols apply to the case of absorbing 
aerosols (m » 1.55-t0.05) as well, except for the addi- 
tional conclusion that the agreements between the 
Rayleigh and Mie results for all models are considerably 
better for the absorbing aerosols than for the nonab- 
sorbing (see Figs. 14-16). The latter conclusion is again 
in line with the one made for monodispersions in Part 
1 . 

D. Error Contours Diagram in — d plane 

Figure 17 shows the error contours for 10% (dashed 
line) and 5% error (solid line) for m ■ 1.55 and models 
M, H and C . The curves indicate an inflection point in 
Haze M (10%) curve and Haze H (5%) at about 9 » 1.5" 
and ^ - 6", respectively, in the range 0.1 nm. 

E. Experimental Design Considerations 

Equation (11) shows how the forwardscattering cor- 
rection factor is included in the transmission law to 
obtain the measured (or apparent) optical depth r'. 
Comments about the two aforementioned experimental 
geometries, discussed in Paper 1, also apply to the case 
of polydisperse aerosols, in addition to the following 
comments regarding the effects of the polydispersity on 
the experimental design. 

In Paper 1 (Fig. 13), results were plotted for /? as a 
function of x9 for the two aforementioned experimental 
geometries. For the polydisperse case, however, R{x9) 
has to be further avemged over the volume extinction 
coefficient to yield R, the path-averaged correction 
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Fi|.l7. Error contour* in tht plane form ■ Min. 

for model* Heze M, Haze H, and Haze C. 


factor for polydisparsa aarosola, dafinad by tha equa- 
tion 


( 18 ) 

The symbols cap( *) and bar(~) denote averaging over 
the pathlength and particle size distribution, respec- 
tively. 

As expressed by Eq. (18), it is not possible to express 
f? as a function * kr^B, which would enable one to 
make some sort of comparison with the results for the 
monodisperse case. However, such a comparison be- 
comes possible if we make the further assumption that 
Q * 2, which is reasonable for the case of large particles 
for which the Rayleigh formulation is valid. Then, Elq. 
(18) reduces to 

dz 2 ^n{t)R{kr„Bz) ^ ^ dzz^n(x)dx. 


where 


( 19 ) 



z - r/r«.z, - ri/r,.Z| - r^/r^. (20) 

ff, obtained by using Eq. (19), can easily be plotted as 
a function of y(* kr^B) for the three size distributions 
for each of two detector systems, as shown in Fig. 18. 
The resulting plots are similar in shape to those for the 
monodisperse aerosols in Paper 1, with the values o(R 
converging to 1 and 0.5 for y -• 0 and ® , respectively. 
However, a comparison between the two sets of plots 
easily shows that increasing the polydispersity results 
in a translation of the curves toward the lower 3' values, 
or, in other words, maximum gradients for R occur at 
lower values of y. 


Fig. 18. Plot* of the path-averaged correction factor v* y, where 
y • kr^a/f and kr^R-i/L for the len»-pinhole detector and open de- 
ft ctor *ysuma. respectively, for models Haze M, Haze H, and Haze 
C. X > 0.5& Mm. 



Fig. 19. Error contours in the r - plane for 
m ■ 1.33, X « 0.55 Mm. for the power law size 
distribution. 



Fig. 20. Error contours in the » - 9 plane for 
m » 1.55-1(0.05), X ■ 0.56 urn. for the power 
law size distribution. 



Fig. 21. Plots of correction factor K vt t) for 
m ■ 1.55-<(0.05), X » 0.55 Mm. and five values 
of the power » for a power law size distribu- 
tion. (Symbols without crosses represent 
Mie results: symbols with crosses, approxi- 
mate resulu.) 
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In order to make acciirate transmission measure- 
ments, it is important that the experimental design be 
based on those values of y for which either R -* lor R 
-* 0.5. But to be able to do so, some reasonable prior 
knowledge of both the mode radius r„ and the polydi* 
spersity is required. 

P. A-Oependence of R and E 

A.1 the results shown in the various plots were com- 
puted for X » 0.55 nm. For any other value X', these 
same residts are valid for the three models. Haze M, 
Haze H, and Haze C, provided r„ is replaced by r'„ such 
that 

rjy - rj\. (21) 

A more detailed parametric study of the behavior of 
the forwardscattering effects involved in different 
transmission measurement experiments designed to 
measure optical extinction and visibility in the atmo- 
sphere will be presented in a separate publication. 

Q. Results for Power Law Size Distributions 

In Fig. 19-21, we present a few selected results for the 
case of a power law size distribution. Figures 19 and 20 
show a series of error contours in ^ plane for refrac- 
tive indices of 1.33-i0.0 and 1.55-i0.05, respectively. 
These two refractive indices usually produced the 
largest and smallest errors (respectively) for a given d 
- V combination. Since an increase in v leads to a re- 
duction in large particle content, the shape of these plots 
is inverted compared with those in Fig. 17. 

Figure 21 shows plots of E vs 0 for a refractive index 
of 1.55-i0.05 and for five selected values of ¥. The 
Rayleigh approximation results are included for com- 
parison. (Again, symbols without crosses represent Mie 
results; those with crosses represent Rayleigh results.) 
The shape of these curves is similar to those in Figs. 
11-16. We see from Fig. 21 that the Rayleigh approx- 
imation is good for this complex refractive index over 
the full range of » values considered. As in the case of 


the Deirmendjian models, the Rayleigh approximation 
is not as good for a real refractive index. 

VIII. Summary Remarks 

From the foregoing discussion, it is obvious that the 
transmission law as expressed by Eq. (1) cannot be used 
except in the case of small particles, for which it is 
merely a good approximation. The effect of not only 
the particle size distribution, refractive index, and 
shape, but also the geometry of the optical system on the 
apparent volume scattering and extinction coefficients 
must be considered. 

This study shows that the R(E) decreases (increases) 
as the size distribution becomes more polydisperse. 
This result can be explained by the fact that as the po- 
lydispersity becomes greater, the relative proportion of 
larger particles increases, resulting in an increase in 
forwardscattered radiation within narrower angles. 
The results are valid for conditions in which multiple - 
scattering effects can be neglected, and only single 
scattering predominates. 

The lens-pinhole detector geometry '/ields the most 
accurate transmission/extinction results, provided the 
desi^ conforms to the prescription that the value of y 
in Fig. 18 should be such that R is close to either 1.0 or 

0. 5. 

This work was supported by NASA grant NSG-1252 
and contract NASl-15198. 
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Table 1. Summary table for size distribution models. 


I 

f»M«r Lw 


)M«I I 

t»Hnl«rU«4 tmMt Uw 


HmI*I > 


■U> 








P, Pi N 


p » •iip(-p,r I 


1/P, 

'i I 11 ♦ P,*P4-iy 


«P^,P4> 


l/P. 


■U J 


Vl 1*P,«P4-1I 


'j 'j P, 

p, (p,-i) a*p,(p.>iti . f-2-, 

A _1 i_i I p,p. 


p, ' 


r •• 0 , 
«(tl ■ 






P, P,-l 




f < 
■(») - 




p, -<i‘p,(p.-m 

-i Ir/B ) ■’ 


p. 1 -p, 1 -p, 

«<rl — ir (r * -r- ) 

P -i * 


r ^ 0, 
lt(rl • 




p.-i 


p, ( 1 -P 4 I 
P,(l*<r/p,) ’) 


x__ » 

p,(p;-TT 


p,(l-(|* 4 -lMr/p,»^) 

?7«vT> 


'‘it* 

'’3' » ' V* ' 




^3< Pa 


P, f^4 


«(«) - 


P, 1-P, 




«r •) 


P, <^/P,) 


-P,»P4-II 


P,(P4-1) 


“ |(P4> 


ORIGINAL PAQt W 
OF. POOR QUAUni 


(continued) 


47 


Table 1. Summary table for size distribution models 
(continued) . 
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Table 1. Summary table for size distribution models 
(concluded) . 
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Figure 1. Relative intensities of sky radiance 
and Sun's direct radiation; r^ rep- 
resents distance from the Sun in solar 
radii (from ref. 14). 
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Figure 2. Geometry of the single 
scattering problem. 
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Figure 3. Schematic illustration of the 
solar aureole photographic 
equipment. 
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Figure 4. (a) A typical solar aureole 3S-mm photograph 

take through a wavelength filter, X = 50 nm. 

(b) An. isodensity tracing of the photograph. 

(c) The computed solar isophote mapping of 
the photograph. 
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Figure 6. The shapes of the conics for various solar 
zenith angles for lens focal lengths 

of 35 mm and 55 mm. 
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Figure 8. Typical size distribution curves 
for parameter v * 4.0, 4.4, and 

5.2. 



Figure 9. A set of phase function cu^es 

corresponding to the size distri- 
bution curves shown in figure 8 
for V values in the range of 
4.0 to 5.2. 
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(b) 


Figure 10. A typical example of inversion of 

solar aureole data to retrieve aero- 
sol size distribution, (a) Solid 
lines denote the final steps to the 
experimental data denoted by the 
samples, and (b) plots represent 
the corresponding retrieved results. 
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(a) 



(c) 


Figure 11. Computer-generated SS circiomsolar 
isophotes for different values of 
the size distribution parameter v: 
(a) 4.0, (b) 4.4, and (c) 5.0. 
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Figure 12. 


Error contours in x-9 
plane for nv = 1.50 and 
percentage errors of 1, 
5, and 10 percent. 
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Figure 13. Aerosol extinction coefficient altitude profile. 
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Figure 14. Model 2 (regularized power law) for n(r). 

Parameter set 2.1: p? variable, P 3 = 2.0, 







